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The proportions of species in many of the five butterfly families (Hesperiidae, Papilionidae, Pieridae, 
Lycaenidae and Nymphalidae) found across the Indian subcontinent show a relatively invariant relationship 
with the overall butterfly species richness, at both local and regional scales. This relationship suggests that 
it is possible to use the species total of a single butterfly family best suited to estimate the overall species 
richness of all other butterflies in an area. Family Papilionidae is a logical choice over others for ease of 
sampling. Also, there is a strong positive correlation between Papilionidae species richness and the overall 
species richness of all other butterflies across all other areas, and the proportion of this family is reasonably 
invariant. The mean proportion (7%) of this family can thus be used to estimate the overall butterfly species 
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richness of an area across the Indian subcontinent for which the Papilionidae species total is known. 
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INTRODUCTION 


The Indian subcontinent, which includes the area from 
Baluchistan (Pakistan) eastwards through India up to 
Myanmar and Sri Lanka, as well as the higher trans- 
Himalayan zone, is habitat for more than 1,439 species 
(Evans 1932, Haribal 1992) of butterflies representing 7.2- 
11.1% of the total world species [13,000 (Owen 1971) - 
20,000 (Vane-Wright 1978)]. Amongst these, about 100 
species are endemic to the Subcontinent (Smetacek 1996) 
and at least 26 taxa are today “globally threatened” (IUCN 
1990). Identification and prioritisation of areas of 
conservation concern, i.e. butterfly biodiversity hotspots, 
are usually based on local endemic and relict taxa, their 
biogeographical affinities and globally threatened and rare 
status. However, prioritisation and selection of such areas 
requires estimation of various ecological indices [e.g. 
Shannon diversity index, Pielou’s evenness index, Similarity 
index (Ludwig and Reynolds 1988)] which depend on the 
‘absolute species richness’ of the species of the area. Data 
on the absolute butterfly species richness of most areas 
across the Subcontinent is non-existent. Traditional methods 
of deriving species richness by collecting and counting all 
the species in an area require much time, effort and 
resources, which were not easily available, and hence 
such studies have not been carried out in India. There is a 
need to evolve easy and cost effective methods to estimate 
the butterfly species richness of areas of concern. 

Beccaloni and Gaston (1995) have proposed such 
a method to predict butterfly species richness of areas 


in the lesser known tropical forests of Central and South 
America, with the help of known species totals of only 
a single sub-family (Nymphalidae: [thomtinae) also 
called the indicator group. This method is based on 
the fact that the proportions of species in many of the 
14 subfamilies and families occurring in these forests 
show a relatively invariant relationship with the overall 
species richness of the area, on both local as well as 
regional scales. Besides, the species richness of this 
indicator group also has a strong positive correlation with 
the overall species richness of all the butterflies across 
the areas and the proportion (4.596) of this group is 
reasonably invariant across tropical forests of central and 
south America. Keeping in mind the findings of Beccaloni 
and Gaston (1995), the present study was conducted to 
determine if proportions of butterfly species in families 
distributed over the Indian subcontinent are also invariant 
with respect to (1) species richness, (11) spatial scale, 
(iii) forest type and (iv) butterfly subregional distribution 
in the Subcontinent. This study further tries to determine 
the potential indicator group amongst the 5 major families 
found in the subregion that can be used to estimate the 
species richness of other butterflies found in different 
areas in the Subcontinent. 


METHODS 
The Indian subregions 


The Indian subcontinent (study area) forms a major 
part of the Oriental region, occupying its extreme 
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northwestern limits. It has been divided into 9 butterfly | Lanka or SL, (ix) Andaman & Nicobar Islands or A&N 
subregions (Evans 1932, Wynter-Blyth 1957) as in of India. Butterfly species found in Baluchistan, Chitral 
Fig. 1. (i) Baluchistan or BA (northern limit up to Safed and higher reaches (1,000-5,100 m) of the Himalaya (east, 
Koh: 26° 00'-34° 00' N and 62° 00'-70° 00' E) and (ii) central and west) have strong Palaearctic affinities 
Chitral or CL (72° 00' E and 36° 00' N) including Chitral, | (Central Asian and Chinese subdivisions), whereas 
Hunza, Baltistan and Ladakh in both Pakistan and India ^ butterflies found in the Peninsular Indian, Malaysian and 
(iii) western Himalayas or WH (Kashmir: 74? 00' Eand — Indo-Chinese subdivisions have strong Oriental affinities. 
36° 00' N to Kumaon: 80° 00' E and 29° 00' N) in India, The drier low-lying areas of PI (Deccan and 
(iv) Central Himalayas or CH (Nepal: 80° 00'-88° 00'E Indogangetic plains) also show affinity with the African 
to 30° 00'-27° 00' N), (v) Northeast India (includes region (Evans 1932, Wynter-Blyth 1957). 

eastern Himalayas from Sikkim: 88° 00' E and 29° 00' N 

to Arunachal Pradesh: 96° 00' E and 28° 00' N in India Methodology 

through Bhutan and parts of Bangladesh) and north Species totals of all the 5 commonly recognised 
Myanmar (up to Shan States: 97° 00'-100° 00'E to butterfly families [Hesperiidae, Papilionidae, Pieridae, 
28° 00'-20° 00' N) or NEI & NM, (vi) south Myanmar Lycaenidae and Nymphalidae (Ackery 1984)]. found 
or SM (Karenni Hills: 97° 00' E and 19° 00' N to Victoria in the Indian subcontinent were gathered. Literature 
Point: 98° 00' E and 10° 00' N), (vii) Peninsular India spanning 117 years (1882-1999) across 69 areas of the 
or PI (Plains and hills of India south of the Himalayas, Subcontinent was reviewed. However, in this paper, 
east of the Indus and west of Brahmaputra), (viii) Sri familywise species totals of only 56 areas, collected 
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Fig. 1: Nine butterfly subregions of the Indian subcontinent and locations of collection sites 
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from 45 different sources of literature, have been used, 
as only these were based on comprehensive surveys 
for which (1) all the 5 butterfly families were sampled 
(i1) had a minimum collecting effort of 2 2 years (51 
areas), (iii) showed no preference to a particular group 
for collection, (iv) covered all the 9 Indian butterfly 
subregions and (v) all the 14 major forest types 
(Champion and Seth 1968) found in the Subcontinent. 
The scientific names of butterflies used in the old records 
were updated and the species correctly placed in their 
respective families, based on the new nomenclature 
(Ackery 1984). The data was then analysed to derive 
the proportions of butterfly species in the 5 families from 
the 56 sites. The areas from which the data were used 
varied from smallest to the biggest site (sites < districts 
< states < sub-regions < the entire Subcontinent), and 
have been ranked on a spatial scale of 1-7, in an 
increasing order (Table 1). The details ofthe areas, their 
relative size, major vegetation types, collecting 
(sampling) effort and source of information are given 
in Table 1. Data on the number of butterfly species per 
family found in each of these areas is summarised in 
Table 2. [For one site “Khasia and Jaintia hills” in 
northeast India, the species totals of 4 families had been 
published by the authors, leaving out the total for 
Hesperiidae, although collections for all the 5 major 
families were done. Hence, the regional proportion of 
Hesperiidae (22.2%) for northeast India was taken as 
an approximate estimate for this site and added to the 
actual species total (464) of the other 4 families 
collected (which thus represented 77.8% of the total 
butterfly species found in this area) to derive the total 
species richness of this area i.e. 596 species. In this 
study, the smallest area in the Subcontinent was New 
Forest, Dehra Dun (4.40 sq. km), which lies in the 
Tropical Moist Deciduous forest zone of the western 
Himalaya. 

However, the type of data used in this study is 
prone to error, including unequal sampling effort across 
areas. Under-recording of species is likely to affect 
the butterfly totals of the least rich areas more than 
those of the richest areas. At the site level, however, 
under-recording is likely to be greatest at the richest area. 


RESULTS AND DISCUSSION 


Patterns in Species Richness 

The proportions of at least 3 out of the 5 families 
(Papilionidae, Lycaenidae and Nymphalidae) are more 
or less independent of the total species richness, or size 
of the area, or forest type, or butterfly subregion, in the 
Indian subcontinent (Table 2; Figs. 2, 3). The 
comparatively invariant relationship exhibited by these 
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families is simpler than the relationships shown by 
families Hesperiidae and Pieridae. The proportion of 
Hespertidae increases with the total species richness, 
whereas that of Pieridae decreases (Fig. 2). This 
variation for these two families is more pronounced 
across the continuous mountainous subregions 
[Baluchistan-Chitral-Himalayas (western-central- 
eastern)-Hills of Myanmar (north-south)] (Fig 1). The 
proportion of Hesperiidae, in general, increases from 
Baluchistan towards south Myanmar [BA(11.8) - CL 
(8.4) - WH (15.1) -— CH (20.1) - NEI & NM (21.9) - 
SM (24.0)], whereas that of Pieridae decreases across 
the same region [BA (21.8) CL (18.1) - WH (10.1) - 
CH (9.9) - NEI & NM (5.9) - SM (5.1); Table 2]. 

As proportions of the first 3 families are invariant 
and show a simple relation to the total species richness, 
it is possible to use the known species totals of the most 
suitable of these three groups in an area to estimate the 
total butterfly species richness of that area. Also, none 
of these 3 groups show ‘saturation’ (Beccaloni and 
Gaston 1995), as their proportions do not decrease with 
the increase in total species richness. Therefore, all three 
are potential indicator groups. 


Selecting an indicator group 

For a group to be an indicator, there should be low 
variance in the relationship between the species richness 
of this group and that of the group we wish to predict 
(Beccaloni and Gaston 1995). Amongst the 3 families 
identified as potential indicators, Papilionidae (x ^ 7.0301; 
SD = 1.1879; n = 51; CV = 16.90) and Lycaenidae 
(x =29.0151; SD = 3.6779; n = 51; CV = 12.68) have 
low variance values (s? < X) for proportion (arcsin 
transformed) of species in families across the 
Subcontinent, as compared to Nymphalidae (X = 33.4740; 
SD = 05.9583; n = 51; CV = 17.80) which exhibits a 
comparatively large variance (s? > X) across the same 
region. Thus, families Lycaenidae and Papilionidae are 
more suitable potential indicator groups than 
Nymphalidae for predicting species richness across the 
Subcontinent. 


Why choose Papilionidae over Lycaenidae as 
indicators? 

Papilionidae (commonly called Swallowtails) are 
taxonomically and ecologically well known in the Indian 
subcontinent, and the distribution of practically all the 
species is known. In contrast, many of the species in 
Family Lycaenidae are very difficult to identify and very 
little is known about their life history and ecology 
because of their obscure habits. Swallowtails (as the 
name suggests, most of them have tails on their hind 
wings) are (i) large in size (wing span: 5 - 19 cm for 
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Fig. 2: Relationship between the proportions of species in a given butterfly family (arcsine transformed) and 
overall butterfly species richness of collection sites (n=56) across the Indian subcontinent 


Indian species) (i1) very active and strong fliers during 
daylight, when they can be observed flying, or feeding 
on flower nectar, or mud puddling, but seldom concealing 
themselves in foliage or settling down to rest, and are 
also (111) eyecatching and colourful, with contrasting 
black as their base colour. In contrast, most of the 
Lycaenidae are (1) very small or medium sized (wing 
span: 1.5 - 6 cm. for Indian species) (ii) cannot be easily 
identified in flight or even at rest, as allied species ofthe 
same genus have similar patterns on the underside of 
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the wings and (iii) are not active fliers like Papilionidae, 
as they are unable to fly for long stretches and soon 
settle down to rest (Haribal 1992, Wynter-Blyth 1957). 
Also, Papilionidae with 94 species (7.01%) is also a 
smaller group to monitor than Lycaenidae with 459 
species (29.22%) (Table 2). All these unique characters 
of Papilionidae make it an easier group to observe, 
identify and sample than Lycaenidae. 

Besides, the average life span of adult Indian 
Papilionids ranges from 20-30 days to a maximum of 
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9 Hesperidae(r = 0.32, p=0.15) 
zu —322. Papilionidae (r * - 0.09, pz0.53) 
— - — - Pieridae (r = - 0.35, p=0.01) 
Lycaenidae (r = 0.15, p=0.28) 
= = = Nymphalidae (r = 0.09, p=0.51) 
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Fig. 3: Relationship between the proportion (%) of species 
in different families of butterflies and the spatial distribution 
of collection sites (n=56). Geographical scale (sq. km) 
1. <10, 2. >10 to 100, 3. >100 to 1000, 4. >1000 to 10,000, 
5. >10,000 to 1,00,000, 6. >1,00,000 to 10,00,000, 
7. >10,00,000 


4 months (Haribal 1992). Their flight period in the plains 
ranges from January to December with many 
overlapping generations, whereas in the hills they fly 
during summer, between April to September, and have 
1-3 generations (Wynter-Blyth 1957), thus Papilionidae 
can be sampled for a longer period in the year. 
Papilionidae are found in all types of habitats (gardens, 
forests, open areas, etc.) from the low lying Indian plains 
to as high as 5,100 m above msl in the Himalaya 
(Common Blue Apollo Parnassius hardwickei Gray) 


Number of other butterfly species 


Number of Papilionidae species 


Fig. 4: Plot of butterfly species richness (excluding 
Papilionidae) versus Papilionidae species richness for sites 
(n=56) across the Indian subcontinent 
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Fig. 5: Distribution of the Papilionidae proportions of different 
collection sites (n=56) across the Indian subcontinent 


(Wynter-Blyth 1957). A large proportion (14.5%) ofthe 
worldwide total of 650 Papilionid species is known to 
occur in the Indian subcontinent (Haribal 1992). 
Papilionid species richness in the Indian subcontinent 
peaks in NEI & NM where a large concentration is 
found [Sikkim (55 species in 7,299 sq. km: Haribal 1992); 
North-east India (62 species in 3,68,000 sq. km) and 
Myanmar (66 species in 6,76,577 sq. km): Wynter-Blyth 
1957}; 

Are Papilionidae good indicators for predicting 
species richness? 

A strong positive relationship exists between 
Papilionidae species richness and the overall species 
richness (of all the other butterfly families) across 
56 different areas over the entire Indian subcontinent, 
and varying on different spatial scales [(1-7; Table 1) 
and (r = 0.980, n —56, p « 0.01; Fig. 4). The histogram 


100 
90 
80 
70 
60 


50 


Expected 


40 
30 
20 
10 





Observed 


Fig. 6: Relationship between the number of Papilionidae 
observed and expected to occur in different sites across the 
Subcontinent (expected values based on the assumption 
that Papilionidae constitute an invariant proportion of the 
total butterfly species found across the entire Subcontinent) 
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Table 2: Total butterfly species richness of areas in the Indian subcontinent and the proportion of species recorded in the families 


Sl.no. Area Total Percentage of total 
Species ‘Hesperiidae Papilionidae Pieridae Lycaenidae Nymphalidae 
ak Indian subcontinent 1439 21.3 (307) 06.5 (94) 06.9 (99) 31.8 (459) 33.3 (480) 
2 Baluchistan 119 11.8 (14) 05.0 (6) 21.8 (26) 28.6 (34) 32.8 (39) 
3 Chitral 166 08.4 (14) 06.6 (11) 18.1 (30) 27.7 (46) 39.2 (65) 
4. Chitral 139 10.8 04.3 15.8 28.8 40.3 
5 Western Himalaya 417 15.1 (63) 07.4 (31) 10.1(42) — 30.9 (129) 36.5 (152) 
6. Kangra Hills 228 110 10.1 16.2 24.6 38.2 
7. Shimla Hills 299 137 07.0 11.4 29.4 38.5 
8. Dehra Dun Valley (New Forest) 148 14.9 07.4 12.8 28.4 36.5 
9. Mussoorne Hills and adjoining areas 323 16.7 57.1 09.9 21.9 39.0 
10.  Mussoorrie Town 146 09.6 06.8 08.9 30.1 44.5 
11. Kumaon Hills 371 14.0 07.0 10.0 29.4 39.6 
12. Central Himalayas (Nepal) 623 20.1 (125) 06.9 (43) 07.9(49) 29.7 (186) 35.3 (220) 
13.  North-East India * North Myanmar 962 21.9 (211) 07.2 (69) 05.9(57) 29.5 (283) 35.5 (342) 
14. NortheastIndia 853 22.2 07.3 06.1 30.2 34.2 
15. Gin 690 230 08.0 07.4 23.5 38.1 
16. Darjeeling distnct 262 10.3 11.1 12.2 18.3 48.1 
17. Naga Hilis 423 15.9 09.0 07.1 26.0 41.8 
18. | Manipur and Naga Hills 321 3754 05.6 00.3 392 17.8 
19.  KhasiaHills 510 192 08.2 07.0 25.0 40.5 
20. Khasia and Jaintia Hills 596 22.2 08.2 06.7 27.8 34.7 
21. = Chin- Lushai (Mizoram) Hills 276 21.4 04.7 09.4 27.9 95.5 
22  N.Chin and Upper Chindwin district — 320 20.3 07.5 08.1 253 38.8 
23. Arakan Coast 159 132 09.4 13.2 26.4 SE 
24. Shan States 228 12:3 07.0 11.0 28.1 41.7 
25. South Myanmar 788 24.0 (189) 06.3 (50) 05.1(40)  34.5(272) 30.1 (237) 
26. Upper Tenassnum 252 19.8 05.6 10.7 26.2 37 7 
27. Tavoy district 401 12.0 08.0 07.0 29:1 37.9 
28.  Merguiandits Archipelago 208 11.5 07.2 13.5 25.5 42.3 
29. Myanmar 1039 25.6 06.3 04.2 32.5 3159 
30.  Peninsular India 315 23.5 (74) 06.0 (19) 10.8 (34) 28.6 (90) 31.1 (98) 
31. Central Provinces 147 21.1 06.1 13.6 23:8 35.4 
32. Calcutta 167 19.2 06.0 11.4 34.1 29.3 
33. South Bihar 124 14.5 07.3 09.7 3/.1 315 
34. | North Bihar 151 129 06.8 15.0 25.2 40. 1 
35. Lucknow 109 138 06.4 17.4 30:3 32.1 
36. Delhi 77 14.3 05.2 26.0 28.6 26.0 
37. Mhow 110 10.9 03.6 26.4 28.2 30.9 
38. Jodhpur ard Mount Abu 78 17.9 07.7 26.9 25.6 21.8 
39. Karachi 70 12.9 04.3 41.4 18.6 22.9 
40. Sind Province 59 16.9 05.1 23.7 32.2 22.1 
41. Central Gujarat (Kathiawar) 78 14.1 06.4 24.4 28.2 26.9 
42. North Gujarat (Kaira) 59 06.8 08.5 309 28.8 25.4 
43. South Gujarat 145 15.2 06.9 17.2 Ska 29.0 
44.  Bombay-Deccan (Pune) 164 11.6 06.1 30.5 25.0 26.8 
45. Konkan Coast 130 a 08.5 14.6 26.9 32.3 
46. N. Kanara district 233 24.0 07.3 09.4 29.6 29.6 
47. Coorg 278 21.6 06.5 11.2 28.8 29.9 
48. Bangalore district 140 14.0 06.4 17.9 35.7 25.7 
49.  Travancore 220 16.8 06.8 11.8 264 37.3 
50. PalniHills 249 18.5 06.0 12.0 33:3 30.1 
51.  Nlilgiri Hills 294 21.8 07.8 11.9 27.9 30.6 
92. Nagalapuram Hills (Eastern Ghats) 117 19.7 10:3 17.9 26.5 25.6 
53. Secunderabad 70 10.0 07.1 20.0 28.6 34.3 
54. Srilanka 242 19.8 (48) 06.2 (15) 12.0 (29) 33.9(82) 28.1 (68) 
55. Andaman and Nicobar Islands 217 19.4 (42) 06.0 (13) 09.2 (20) 31.3 (68) 34.1 (74) 
56. Great Nicobar Island 68 07.3 08.8 14.7 32.4 36.8 


Values in parenthesis are species totals for nine butterfly sub-regions and the whole of the Indian subcontinent so far known 
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(Fig. 5) demonstrates that the variance of Papilionidae 
proportions across these areas is reasonably low, with 
Papilionidae constituting 6-8% of the butterfly species 
in 28 out of 56 areas (and 6.5-7.5% in 18 out of 
56 areas) in the Subcontinent. 

Five areas [Mhow, Chitral, Chin-Lushai (Mizo) 
Hills, Darjeeling Hills and Nagalapuram (Nagari) Hills; 
Table 2] could have been excluded from this analysis, 
thereby increasing the level of correlation. The first three 
have the lowest proportions of Papilionidae (Mhow 3.6%, 
Chitral 4.3% and Chin Lushai 4.7%) of all areas. Data 
of Mhow and Chin-Lushai suffer from sampling error 
as less than one year of collecting effort was 
undertaken. However, data from Chitral also includes a 
large number of Palaearctic species besides the Oriental 
species (as this region has strong affinities with the 
Palaearctic region), which are likely to decrease the 
proportion of Oriental species. The last two areas, on 
the other hand, have the highest proportions of 
Papilionidae of all areas [Nagalapuram Hills (10.3%) 
and Darjeeling (11.1%)]. Data from Nagalapuram Hills 
(with only 1.5 years of sampling) is also under-sampled, 
particularly for Nymphalidae (Table 2). On the other 
hand, the exact sampling period for Darjeeling district is 
not mentioned in the original text (Maude 1949) but the 
data reflects low sampling of species from this area, 
particularly those of the families Hesperiidae and 
Lycaenidae (Table 1). 

A combined data set for sites, districts, states and 
regions was tested against random draw model in which 
the proportion of Papilionidae in each area was assumed 
to equal that for the whole of the Subcontinent (6.5%). 
Correlation between the number of Papilionid species 


observed and expected to occur in each area is high 
(r = 0.982; n = 55, p= 0.001) and the relationships are 
fitted by a slope of | (Fig. 6). 


CONCLUSION 


A reasonably invariant relationship exists between 
proportions of Papilionid species and the overall butterfly 
species richness across the Indian subcontinent, 
independent of (i) the different butterfly sub-regions 
(11) forest types found in the region (111) different spatial 
scales and (iv) species richness of areas. This suggests 
that it should be possible to use the average proportion 
(7%) of this family to estimate the butterfly species 
richness of areas across the Indian subcontinent, for 
which the Papilionidae species total is known. Using 
the figure for the overall butterfly richness of an area, 
approximate estimates for the species totals of the 
individual butterfly groups in the area can be gene- 
rated on the basis of known relationships between total 
butterfly species richness and the species richness of 
the groups across all areas (provided that these are 
strong). In the Indian subcontinent, the species rich- 
ness of at least three families, Papilionidae, Lycaenidae 
and Nymphalidae, show little variation, and their 
regional proportions can be used to estimate the species 
totals. However, for families such as Hespertidae and 
Pieridae, that show a larger variance across the 
9 butterfly sub-regions in the Indian subcontinent, sub- 
regional proportions are already known (Table 2) and 
would give more precise estimates of this relationship 
for predicting the species richness of the respective 
families. 
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